AD-A158  0$4 


THESIS 


Approved  for  public  release,"  distribution  is  unlimited 


85  "  8  12  065 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (V**"  D«»  Kntmrmd) 


REPORT  DOCUMENTATION  PAGE 


.  REPORT  NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT’S  CATALOG  NUMBER 


4.  TITLE  (and  Subthlm) 

Localization  of  Buried  Objects  in  Water- 
Saturaxed  Sand  by  Variable  Incidence 
Acoustic  Pulse  Reflections 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Master's  Thesis 
June  1985 


«.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORfF) 


Gwan-Sik  Bang 


•  .  CONTRACT  OR  GRANT  NUMBERS 


>■  PERFCRMING  ORGANIZATION  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93943-5100 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 


<1.  CONTROLLING  OFFICE  NAME  AN.  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93943-5100 


I 


MONITORING  AGENCY  N /.ME  4  AOORESSf//  dllloront  tram  Controlling  Ottleo) 


12.  REPORT  DATE 

June  1985 


13.  pages 

56 


IS.  SECURITY  CLASS,  (at  th to  report.) 

Unclassified 


15 m.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  'oi  thi»  Riport) 


Approved  for  public  release;  distribution  is  unlimited 


17.  DISTRIBUTION  STATEMENT  (at  tho  mbotrncl  ontorod  In  Block  30,  II  dllloront  tram  Report ) 


19.  KEY  WORDS  (Continue  on  r,ovr»o  rid*  It  mnd  tdmntlty  by  block  numb  or) 

2  2 

T  -X  Method;  Ray'  Parameter  Method,  Acoustic  Pulse  Reflections; 
Water-Saturated  Sand 


ZO.  ABSTRACT  (Conttmjo  on  reeeree  tldt  It  neeeaeary  and  tdonuty  by  block  namborf 

The  purpose  of  this  project  is  to  localize  buried  objects 
in  water-saturated  sand  by  variable  incidence  acoustic  pulse 
reflections.  In  particular,  thin  aluminum  plates ’.are  used  to 
model  the  buried  objects.  The  sound  speed  and  the  depth  of  the 
overlying  sand  layer  ace  predicted  by  using  the  T2«X2  and  the 
ray  parameter  methods.  The  experimentally  computed  results 


DD  t  1473  EDITION  OF  I  NOV  *1  It  OBSOLETE 

S  N  010  2-  LF-  01  41  6401 


UNCLASSIFIED 


UNCLASSIFIED 


security  classification  of  this  rage  p*t>w«»  om»  shw« 


#20  -  ABSTRACT  -  (CONTINUED) 

have  an  average  3.5%  depth  error.  The  calculated  speed 
of  sound  in  the  water-saturatsd  sand  layer  differed 
by  8.1%  between  the  t2-x2  and  the  ray  parameter 
methods. 


S  N  0102-  LF.0U-6601 


UNCLASSIFIED. 


SECURITY  CLASSIFICATION  OF  THIS  FAflCfOftaa  Oar*  Sni*rw0 


Approved  for  public  release;  distribution  is  unlimited 

Localization  of  Buried  Objects  in  Water-Saturated  Sand 
by  Variable  Incidence  Acoustic  Pulse  Reflections 


Gwan-Sik  Bang 

Major,  Republic  of  Korea  Army 
B.S.,  Han  Yang  University,  1975 


Submitted  in  partial  fulfillment  of  the  |  Aocno-s‘rT 
requirements  for  the  degree  of  r,—  -.  r~ 


MASTER  OF  SCIENCE  IN  PHYSICS 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1985 


By - 

Dir-trVv. 

Avr.il-'-' 


Author : 


Approved  by ; 


^  Gyan-Sik  Bahg 

a 

l  Suk  Wang  Yoortl  Tnfesis  Advisor 

i  ■ 


Thomas  B/ 


-ielson,  Co-Advisor 


Gordon  E.  Schacher,  Chairman, 
Department  of  Physics 

. <d.  ji  , 

7  John  N.  Dyer, 

Dean  of  Science  and  Engineering 


ABSTRACT 


The  purpose  of  this  project  is  to  localize  buried 

objects  in  water-saturated  sand  by  variable  incidence 

acoustic  pulse  reflections.  In  particular,  thin  aluminum 

plates  are  used  to  model  the  buried  objects.  The  sound 

speed  and  the  depth  of  the  overlying  sand  layer  are  predicted 
2  2 

by  using  the^T  -X  and  the  ray  parameter  methods.  The 
experimentally  computed  results  have  an  average  3.5%  depth 
error.  The  calculated  speed  of  sound  in  the  water-saturated 

2  2 ' 

sand  layer  differed  by  8.1%  between  the  T  -X  and  the  ray 
parameter  methods.  r'  ,  j  ;  >  -j  <;  \  :<  ,'r  •<:_  V ■’  '■  ’’  f'”  '  ' 
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I. 


INTRODUCTION 


A.  OBJECTIVES 

Identification  of  buried  objects  in  the  ocean  bottom  is 
an  important  problem  in  many  fields  such  as  oil  exploration, 
geophysical  research,  and  ocean  engineering.  A  simpler 
problem  is  the  identification  of  layers  of  infinite  extent 
and  unknown  properties  in  the  ocean  sediments.  The  objec¬ 
tives  of  this  thesis  ere: 

1.  to  determine  experimentally  the  depth  of  burial 
of  a.  thin  plate  in  a  simulated  sedimentary  bottom; 

2 .  to  determine  the  sound  speed  and  the  thickness 
of  the  overlying  sand  layer;  and, 

3 .  to  determine  the  sound  speed  and  thickness  of 
the  thin  plate. 

It  was  the  intent  of  this  thesis  to  determine  all  of  these 
properties  by  remote  sensing  with  acoustic  pulses  in  order 
to  simulate  an  ocean  measurement  with  a  ship  and  hear- 
surf ace  sensors . 

The  first  objective  could  have  been  accomplished  by 
normal-incidence  time-of-f iight  measurements;  however,  the 
second  objective  requires  more  information  t.ian  that  so  a 
variable  angle  of  incidence  technique  was  used.  This  second 
objective  was  actually  the  principal  goal  since  the  third' 
objective  could  not  be  achieved  by  using  the  technique 
described  herein.  The  technique  failed  because  the  pulse 
could  not  be  made  short  enough  to  distinguish  the  reflection 
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from  the  plate's  lower  surface.  This  is  primarily  the 
result  of  using  very  thin  plates  with  respect  to  the  acoustic 
wavelength.  Perhaps  the  amplitude  of  the  reflection  from 
the  top  of  the  plate  could  be  used  to  assist  in  identifying 
the  plate  material  but  the  overlying  sand  layer  makes  this 
very  difficult. 


B .  APPROACH 

The  approach  used  is  as  follows: 

1.  Short  pulses  consisting  of  several  cycles  of  a 
15CkHz  sinusoidal  wave  were  used  as  the  source 
pulse. 

2.  A  directional  projector  was  used  to  minimize  inter¬ 
ference  from  surface  and  side-wall  reflections. 

3.  The  receiving  hydrophone  was  mounted  on  a  carriage 
so  that  the  horizontal  distance  from  source  to 
receiver  could  be  easily  varied. 

4.  Each  of  the  thin  plates  were  buried  in  the  sand 
bottom  of  a  small  water  tank.  ' 

5.  The  plates  were  aluminum  ranging  in  thickness,  from 
0.1  to  1.8  cm  and  were  buried  at  depths  from  1.8 
to  5.0  cm. 

'  6.  Two  analytical  methods  were  used  for  reducing  the 
data.  One  method  was  based  on  the  relationship 
between  the  square  of, the  travel  time  (T)  of  the 
reflections  and  the  square  of  the  source-to-receiver 
horizontal  range  (X) .  The  other  method  was  based 
on  maintaining  a  constant  ray  parameter  (P)  defined 
as  the  sine  of  the  incident  angle  divided  by  the 
local  sound  speed. 


C.  BACKGROUND 

Reference  1  gives  an  excellent  historical  development 
of  various  seismological  methods  for  determining  properties 
of  sediment  layers.  The, following  is  a  brief  summary  of 


that  account. 


Karcher  and  Haseman  made  the  first  successful  seismic 


reflection  measurements  in  1921.  Techniques  for  acoustic 

determination  of  sediment  properties  have  advanced  greatly 

since  then  and  generally  fall  into  one  of  three  categories: 

normal  incidence,  wide-angle  reflection,  and  refraction 

methods.  Most  of  the  development  has  been  done  by  the 

petroleum  exploration  industry;  however,  the  geophysical 

community  has  used  these  methods  to  study  the  earth’s  crust 

as  early  as  the  late  lS30's.  Dix,  Krey  and  Durbaum 

refined  the  sound  speed  measurements  for  sediment  layers 

below  the  first  layer  by  properly  accounting  for  the  shape 

of  the  arrival  time  as  a  function  of  range  curves.  3y 

1955,  Dix  had  demonstrated  the  equivalence  of  T-X  methc 
2  2 

and  T  -X  method  and,  in  1968,  Le  Pichon  refined  the 
2  2 

T  -X  method  until  it  became  a  standard  analysis  tool. 

A  major  enhancement  was  introduced  by  Bryan  (see  Ref.  2)  in 

2  2 

1974  with  the  ray  parameter  method.  Both  the  T  -X  method 
and  the  ray  parameter  method  were  used  in  the  data  analysis 
for  this  thesis. 


II.  THEORY 


A.  T-X  AND  T2-X2  METHODS 

A  simple  model  of  a  layered  sediment  is  shown  in  Fig. 
2.1.  The  layers  have  thicknesses  H^ ,  H£  and  H^  and  sound 
speeds  V^,  V7  and  .  The  T2-X2  method  [Ref.  1]  can  be 
applied  to  this  situation  as  described  in  this  section. 


Figure  2.1.  Geometry  of  Layers 

For  convenience,  let  the  sound  speed  in  the  uppermost 
layer  (the  water  layer)  be  less  than  that  in  the  first 
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sediment  layer.  The  length  of  the  p  .th  that  reflects  from 
the  top  of  the  sediment  layer  is  given  by. 


X2  +  4  (Hx  -s)2 


(1) 


To  convert  this  to  travel  time  squared,  divide  through 


by  V2: 


T2(X)  =  T2(0)  +  (-i)X2. 

/'  '  V2 


(2) 


where  T(X)  is  the  reflection  travel  time  at  horizontal  range 

X  and  T  ( 0 )  is  the  intercept  of  the  travel  time  curve  at 
2 

X  =  0.  Since  T  (0)  is  a  constant,  the  reflection  curve  is 

a  hyperbola.  The  corresponding  curve  for  the  direct  path 

from  the  source  to  the  receiver  is  a  straight  line  and  this 

hyperbola  is  asymptotic  to  that  direct-path  line. 

2  2 

When  plotted  in  the  T  -X  plane,  however,  these  lines 

are  both  straight.  In  this  way,  the  slope  of  the  reflection 

curves  can  be  easily  determined.  As  can  be  seen  from  the 

2  2 

equation  for  the  reflection  path,  the  slope  in  the  T  -X 

plane  is  -i-  which,  in  this  thesis  work,  is  calculated  by 
V1  2  2 

linear  regression  of  the  T  -X  data.  This  rather  simple  .  . 
but  accurate  method  is  due  to  Green  {Ref.  1]  and  is  asso¬ 
ciated  with  the  wide-angle  seismic  reflectioh  technique. 

The  paths  that  pass  through  one  or  more  sediment  layers 
before  reflection  (see,  for  example,  R^2  in  Fig.  2.1)  do  not 
lead  exactly  to  hyperbolas  but  the  departure  is  net  great 
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and  only  a  small  amount  of  error  is  introduced  by  fitting 

2  2 

this  data  with  straight  lines  in  the  T  -X  plane. 

If  the  layers  were  extremely  thick  or  had  great  differ¬ 
ences  in  sound  speed,  then  the  error  might  become  apprecia¬ 
ble.  The  reciprocal  of  the  slope  of  the  R^2  curve  is  equaj. 
to  the  average  speed  (squared)  through  both  the  water  and 
the  first  sediment  layer: 

V  =  (H1+H2)/(  (H1/V1.)  +  (H2/V2)) 


=  (v1t1+v2t2)/(t1+t2) 


(3) 


Representative  T-X  curves  for  this  case  are  shown  in 

Fig.  2. 2. a.  Curve  D  represents  the  direct  path  and  R^ 

represents  the  reflection  curve  for  reflection  from  the  top 

of  the  sediment.  For  R12'  t^ie  exact  curve  (solid  line)  is 

given  as  veil  as  the  hyperbola  approximation  (dotted  curve) 

2  2 

to  that  curve.  The  equi valient  curves  in  T  -X  space  are 


shown  in  Fig.  2.2.b.  Here 
speed  squared.  The  interce 
vertical  travel  time  through  the  layers. 

The  principal  limitatic 
in  these  figures.  For  any 

tion  from  the  top  of  the  sejdiment,  the  curve  is  not 
exactly  a  straight  line  in 
slightly  below  the  equival^ 
increases. 


[the  inverse  slopes  give  the  sound 
pt  T ( Q )  gives  the  two-way 


2  2 

n  of  the  T  -X  method  is  shown 


reflection  other  than  the  reflec- 


2  2 

T  -X  coordinates  but  it  falls 


nt  straight  line  as  range 
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return  from  the  bottom  with  that  same  parameter  value 
although  they  arrive  at  the  receiver  at  different  ranges. 

Note  that  these  rays  are  parallel  through  the  water  column. 
The  travel  time  and  horizontal  distance  for  these  two  rays 
are,  according  to  Eqs .  (7)  and  (8), 

AT  (P)  = 

AX  (P)  = 

where  the  factor  2  accounts  for  the  round  trip  of  the  ray. 

If  the  sound  speed  is  constant  at  c^  in  the  layer,  these 
equations  can  be  integrated  to  give, 

-1  2  2  ~1/2 

AT  =  2c,  (1  -Pc)  Az  •  (11) 

2  2 

AX  ■  -  2Pc1(l  -P  cj)  Az  (12) 

The  expressions  for  sound  speed  and  layer  thickness  are 
then , 

C1 

Az 
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,1  AX> 1/2 
'P  AT 


1  2  2 
±cxAT(l  ~P*cp 


1/2 


(1  >) 


(14) 


,Zl+<iZ  -1  2  2  '1/2 
2  /  c  A  (1  -Pc)  dz 


(9) 


/1+AL  2  2  "1/2 
2  /  Pc (1  -  Pc}  dz 


(10) 
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Figure  2.4.  T-X  Curve  and  the  Ray  Parameter 

The  important  feature  of  the  ray  parameter  method  is  that 
the  structure  above  the  layer  has  dropped  out  because  of 
the  choice  of  ray  path. 

In  order  to  apply  this  technique  in  practice,  an  initial 
value  of  ray  parameter  was  chosen  based  on  a  geometry  that 
gives  good  clean  reflection  signals.  After  the  T-X  curves 
are  drawn,  this  value  of  ray  parameter,  which  gives  the 
slope  of  the  T-X  curves  at  the  points  of  interest,  was  used 
to  identify  the  two  parallel  return  rays.  The  principal 
difficulty  of  this  method  is  that  the  T-X  curves  are  not 
straight  lines  and,  therefore,  curves  must  be  fit  to  the 


experimental • data.  This  introduces  some  uncertainty  into 
the  process  but,  fortunately,  the  values  for  layer  thickness 
is  not  very  sensitive  to  the  type  of  fit  that  is  done. 


1XPERIMENTAL  PROCEDURES 


A.  ARRANGEMENT  OF  THE  EXPERIMENT 

•  The  measurements  were  made  in  a  tank  which  was  300  cm 
long,  115  cm.  wide  and  90  cm  deep.  The  tank  was  filled,  with 
care-  ully  de-aerated  (fresh)  water  and  had  a  55  cm  deep 
sand  bottom  (also  de-aerated) . 

At  first,  an  ITC  501  projector  was  used  but  the  beam 
pattern  was  fairly  bread  leading  to  excessive  reflections 
from  the' surface  and  the  side-walls  and,  more  importantly, 
this  projector  was  large  enough  so  that  it  would  ring  for  a 
significant  length  of  time  Compared  to  the  desired  pulse 
length.  As  a  result,  the.  output  pulse  was  too  long  regard¬ 
less  of  the  length  of  the  pulse  used  as  input. 

These  problems  led  to  the  selection  of  a  USRF  F-41 
projector,  which  is  a  directional  laboratory  standard  trans¬ 
ducer  used  in  acoustic  reciprocity  calibrations  over  the 
frequency  range  from  15  to  150  kHz.  This  projector  consists 
of  lead  zirconate-titagate  elements • cemented  to  a  high- 
density  kennametal  disks.  The  resulting  array  of  elements 
is  about  3.8  cm  wide  and  5  cm  high. 

The  receiving  hydrophone  used  was  an  omnidirectional 
LC-10  hydrophone  mounted  on  a  carriage  that  was  free  to  roll 
along  tracks  installed  on  the  long  sides  of  the  tank.  In 
this  way  the  horizontal  distance  (X)  could  be  easily  varied. 
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A  block  diagram  of  the  electronics  used  is  shown  in 
Fig.  3.1.  The  digital  oscilloscope  was  used  to  measure  the 
time  differences  between  various  pulses.  In  order  to  reject 
unwanted  noise,  a  filter  was  used  to  pass  only  100  kHz  to 
200  kHz.  ■ 

The  choice  of  pulse  length  was  governed  by  the  criteria 
that  it  be  (1)  long  enough  to  reach  steady-state  conditions 
in  the  projector  and  (2)  short  enough  to  resolve  the  indi¬ 
vidual  reflections  from  the  layers.  The  best  compromise 
turned  out  to  be  four  cycles  of  a  15,0  kHz  sine  wave,  but, 
unfortunately,  this  was  not  short  enough  to  completely 
define  the  thin  plate.  At  150  kHz,  one  cycle  is  6.7  pS 
so  the  four-cycle  pulse  was  about  27  pS  long.  A  typical 
received  signal  using  the  four-cycle  pulse  is  shown  in 
Fig.  3.2  while  Fig.  3.3  shows  the  same  situation  with  an 
eight-cycle  pulse.  The  four-cycle  pulse  gives  much  better 
definition  of  individual  arrivals. 

Before  the  actual  measurements  were  made,  : 
were  performed  in  order  to  insure  that'  side-wa 
reflections  would  not  interfere  with  the  measu: 
depths  of  the  projector  and  hydrophone  were  de 
trial  and  error  not  only  to  minimize  interferi; 
but  also  to  optimize  the  separation  between  th 
arrivals.  In  addition,  this  set  a  practical  1 
horizontal  range. 

An  attempt  was  made  to  move  the  axis  of  th 
receiver  line  toward  the  tank  wall  so  that  two 


;everal  tests 
.1  and  surface 
remen  ts.  The 
rermined  by 
ig  reflections 
j  desired 
Lmit  on  the 

2  source- to- 
parallel 
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Figvr-a  3.2.  4  Cycle  Pulse  Signal 


Figure  3.3.  8  Cycle  Pulse  Signal 


measurements  could  be  made  without  disturbing  the  sand  to 
bury  each  plate.  This  introduced  unacceptable  side-wall 
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reflections  so,  instead,  three  plates  were  buried  in  line 
along  the  long  axis  of  the  tank  with  the  projector  and 
receiver  position  as  shown  in  Fig.  3.4. 

Here,  the  plates  were  buried  at  different  depths  and  the 
projector  and  hydrophone  were  positioned  as  shown.  (Notice 
that  in  the  third  set-up  on  the  right  rhe  source  and  receiver 
are  reversed  to  avoid  reflections  from  the  end-wall.) 

Figure  3.5  shows  the  side-wall  reflections  that  are 
introduced  when  the  source-to-receiver  axis  is  displaced 
only  15  cm  toward  one  side  of  the  tank.  This  result  should 
be  comparable  to  Fig.  3.2  but  only  a  position  along  the 
centerline  of  the  tank  turned  out  to  be  acceptable. 

B.  EXPERIMENTAL  PROCEDURE 

Once  the  experiment  was  set  up,  the  first  step  was  to 
identify  the  individual  arrivals  from  a  particular  buried 
plate.  By  considering  the  projector  depth  S,  the  hydrophone 
depth  h,  the  water  depth  H,  the  major  incident  angle,  and 
the  slant  distance  R,  the  time  of  arrival  of, the  various 
reflections  was  determined  by  simple  image  theory  [Ref.  3,. 
pp.  427-429] .  A  representative  diagram  of  the  various 
arrivals  is  shown  in  Fig.  3.6  for  constant  sound  speed, in 
the  water  layer. 

Since  a  directional  projector  was  used  and  was  aimed 
downward  toward  the  sediment,  the  surface-reflected  path 
7  could  be  neglected.  For  most  of  the  measurements,  paths  1, 
3  and  4  were  used. 
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PLATE 


Figure  3.6.  Waves  in  the  Isospeed  Shallow. 

Water  Channel 

Basically,  three  methods  were  used  to  initially  identify 
the  arrivals  on  the  digital  oscilloscope.  First,'  the 
arrival  times  predicted  by  the  method  of  images  were  com¬ 
pared  to  the  measured  arrival  times..  The  amplitude  of 
the  arrivals  also  helps  in  the  identification  process,  since 
the  pressure  amplitude  is  reduced  by  the  reflection  coeffi¬ 
cient  at  each  of  the  interfaces  of  the  layers'.  Finally,  a 
path  that  reflects  from  the  surface  can  be  readily  identi¬ 
fied  by  disturbing  the  water's  surface  with  a  ruleri  Any 
surface-reflected  arrivals  change  greatly  on  the  oscillo¬ 
scope  display  and  can  easily  be  separated  from  those 
arrivals  that  'only  reflect  from  the  bottom. 


A  typical  display  on  the  oscilloscope  is  shown  in  Fig. 
3.7.  Here,  the  projector  is  inclined  so  that  the  beam  has 
an  incident  angle  of  55  degrees  and  the  water  depth  is  35  cm. 
The  arrivals  shown  are: 

1.  direct  wave, 

2.  surface-reflected  wave, 

3.  bottom-reflected  wave, 

4.  wave  reflected  from  the  top  of  the  plate, 

5.  wave,  reflected  from  the  bottom  of  the  plate  after 
passing  through  the  plate, 

6.  bottom-surface-reflected  wave. 


1  2  3  4  5  6 


Figure  3.7.  Different  Kinds  of  Waves  on  Oscilloscope 


Once  the  arrivals  are  identified  for  one  source-to-receiver 


range,  the  signals  can  be  followed  by  eye  if  the  receiver 
carriage  is  moved  smoothly. 

Next,  the  range  X  is  varied  by  moving  the  carriage  and 

the  arrival  times  T  of  the  individual  reflections  are 

recorded.  These  measurements  form  the  basis  of  the  T-X  and 
2  2 

T  -X  methods.  When  these  measurements  are  plotted  as 

arrival  time  T  as  a  function  of  range  X,  the  direct  arrivals 

appear  as  straight  lines  and  the  reflected  arrivals  (either 

from  the  top  of  the  sand  layer  or  from  any  of  the  layers 

below)  appear  as  parabolas .  (They  are  actually  closer  to 

hyperbolas  but,  over  the  limited  range  available  in  this 

experiment,  parabolas  were  used  as  an  approximation.) 

2 

The  T  and  X  measurements  were  also  plotted  as  T  as  a 

2  . 
function  of  X  .  On  these  coordinates,  any  of  the  lines  are 

straight  (or  nearly  so) .  From  the  slope  of  these  lines, 

the  sound  speed  either  in  the  water  or  in  the  layers  could 

be  calculated.  A  linear  regression  method  was  used  to 

compute  the  slope  for  the  individual  arrivals. 

Next,  the  ray  parameter  P  was  calculated.  This  value 

is  fixed  for  a  given  measurement.  The  incident  angle  was 

determined,  from  the  source-receiver  geometry  and  the  T-X. 

data  were  fit  with  parabolas.  Then  P  could  be  calculated 

2 

directly  from  the  parabola  equation:  T  =  aX  +  bX  +  c , 
dT/dX  =  P  =  2aX  +  b. 

Once  the  parabolas  have  been  determined,  the  layer  depth 
•*  found  from  the  T-X  measurements.  For  example,  the 


arrival  from  the  wave  that  passed  through  the  upper  sand 

layer  would  be  described  by  a  parabola  of  the  form 
2 

=  a^X^  +  b^X^  +  c^  and  the  reflection  from  the  top  of  the 

2 

sand  layer  would  be  given  by  T2  ~  a2X2  +  b2X2  +  c2*  Tbe 
ray  parameter  P  must  be  the  same  for  each  arrival;  there¬ 
fore,  for  the  wave  passing  through  the  sand  layer, 

P  =  2a^X1  +  b^  or  X.^  =  (P-b^/^a^,  and,  for  the  wave  reflected 
from  the  top  of  the  sand  layer,  P  =  2a2X2  +  b2  or  X2  =  (P-b2)/2a2 
Since  AX  =  X-^-X.,, 


P-b  -  P-b1  P-b„  -  P-b_ 

AT  =  +  bl(_2a^)  +  cl}  "  {a2(_2aJ)  +b2  {~2a^)  +C2 }  ] 


These  two  values  (T  and  X)  are  necessary  to  calculate  the 

2  2 

depth  and  the  sound  speed  (to  compare  with  the  T  -X  deter¬ 
mination  of  sound  speed)  of  the  layer.  The  values  of  layer 
depth  Az  and  sound  speed  c^  are  given  by, 

1  22  1/2 
•Az  =  j  c1AT(l  -P  cp 

2  2 

The  first  value  of  c^  is  from  the  T  -X  method  and  AT  is 
from  the  parabolic  fit  of  the  T-X  data.  The  new  value -of 
c^  can  be  compared  to  the  old  value  as  a  check  on  the  two 
methods . 

C.  SUMMARY  OF  EXPERIMENTAL  PROCEDURE 

1.  Bury  the  thin  homogeneous  plates  at  various  depths? 

2.  Identify  the  arrivals  on  the  oscilloscope; 
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3.  Measure  the  horizontal  distance  X  and  the  arrival 
time  T  of  the  various  arrivals; 

4.  Make  the  T-X  and  T^-X^  plots; 

5.  Find  the  sound  spaed  c ,  in  the  layer  by  using  the 

T5-x2  method;  1 

6.  Determine  the  incident  angle  from  the  source- 
receiver  geometry; 

7.  Calculate  the  ray  parameter  P  (fixed  for  each 
experiment) ; 

8.  Fit  the  T-X  data  with  parabolas; 

9.  Calculate  AX  and  AT  by  the  ray  parameter  method; 

10.  Calculate  the  layer  depth  Az; 

11.  Calculate  c^  and  compare  to  the  value  determined 
previously. 


r~ 
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IV. 


ANALYSIS  AND  RESULTS  OF  EXPERIMENTS 


A.  MEASUREMENT  AND  ANALYSIS 

A  measurement  of  the  sound  speed  in  the  water  was 

required  for  each  of  the  reflection  analyses  so  this  was 

2  2 

measured  using  the  T  -X  method.  The  slope  of  the  straight 
line  fit  to. the  data  points  gave  a  sound  speed  of  1447  m/s 
and  the  correlation  coefficient  associated  with  the  linear 
regression  was  almost  100%.  For  this  measurement,  the  pro¬ 
jector  beam  was  aimed  horizontally  and  the  hydrophone  was 
positioned  at  the  same  depth  as  the  projector. 

In  each  of  the  subsequent  experiments,  the  projector  and 
hydrophone  depths  were  chosen  to  optimize  the  separation 
and  clarity  of  the  reflection  signals.  In  every  case,  the 
source  signal  was  four  cycles  of  a  150  kHz  sine  wave. 

1 .  Measurement  1 

For  this  experiment  an  aluminum  plate  (thickness  of 
i.7  cm)  was  buried  3.5  ±  0.05  cm  below  the  sand's  surface 
and  the  projector  depth  was  8.5  cm  while  the  hydrophone 
depth  wag  6.5  cm. 


TABLE  ,1 

Speed  of  Sound  (Measurement  1) 


WAVE 

SPEED  OF  SOUND  (m/s) 

CORRELATION  (%) 

ERROR  {%) 

SAND 

1661 

98.3 

BOTTOM 

1415 

99*4 

2.2 

DIRECT 

,  1463 

99.8 

1.1 

Table  I  gives  the  speed  of  sound  calculated  from  the  various 
arrival?  along  with  the  correlation  coefficients  of  the 
linear  regressions.  The  speed  of  sound  in  the  water  should 
be  compared  to  the  value  of  1447  m/s  discussed  above.  The 

~  2 

T-X  and  'i  -X  data  are  shown  in  Fig.  4.1. 

As  discussed  in  the  previous  chapter,  the  T-X  data 
were  fit  with  a, parabola  for  simplicity.  The  data  for  the 
wave  that  reflected  from  the  top  of  the  aluminum  plate  was 
fit  with 

T  =  0.0926  X2  -  2.89  X  +  459 

and  the  data  for  the  wave  that  reflected  directly  from  the 
top  of  the  sand  was  fit  with 

T  =  0.113  X2  -  3.03  X  +.  410 

The  geometry  selected  for  the  ray  parameter  determination 
gave  an  incident  angle  cf  29  degrees  which  corresponds  to  a 
ray  parameter  P  equal  to  3.32  xlO  ®  s/m.  This  value  of  slope 
locates  a  point  on  each  of  the  two  parabolas.  The  differ¬ 
ence  in  X  values  gives  a  AX  of  5.47  cm.  The  parabolic  form 
for  the  reflection  from  the  top  of  the  plate  gives  the  time 
difference  AT  of  51.4  yS. 

Next,  the  layer  thickness  Az  was  calculated  from 
Equation  (14).  In  this  case,  Az  was  determined  to  be  3.56  cm 
This  differs  by  only  0.‘3%  from  the  measured  value  for  the 
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r  X  PLOT  OF  MEASUREMENT  1  -  xVijot  of  measurement  i 


overlying  sand  layer.  (The  burial  depth  of  the  plate  gives 
the  thickness  of  the  top  sand  layer.)  The  value  for  the 
layer's  sound  speed  was  also  calculated  from  the  ray  param¬ 
eter  method  but  its  value  was  not  significantly  different 

2  2 

from  that  calculated  by  the  T  -X  method. 

2 .  Measurement  2 

For  the  second  measurement,  the  hydrophone  depth  was 

6.4  cm,  the  projector  depth  was  8.5  cm  and  the  aluminum  plate 

(1.7  cm  thick)  was  buried  at  a  depth  of  2.9  ±  0.1  cm  in  the 

sand.  The  calculated  sound  speeds  are  given  in  Table  II  and 
2  2 

the  T-X  and  T  -X  data  are  shown  in  Fig.  4.2. 


TABLE  II 

Speed  of  Sound  (Measurement  2) 


WAVE 

SPEED  OF  SOUND  vm/s) 

CORRELATION  (%) 

ERROR 

SAND 

.1647 

94.3 

■■ 

BOTTOM 

1419, 

99.6  ' 

1.96 

DIRECT 

1437 

99.1 

0.69 

Figures  4.3  and  4.4  show  representative  displays  from 
the  oscilloscope  at  two  of  the  ranges  at  which  measurements 
were  taken. 

The  direct  wave  (A) ,  the  surface-reflected  wave  (B) , 
the  wave  that  reflects  from  the  top  of  the  sand  (C) ,  the 
wave  that  reflects  from  the  top  of  the  buried  plate  (D) ,  and 
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Figure  4.2.  T'-X  and  T  -X  Plots  (Measurement  2) 
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the  wave  that  reflects  from  the  top  of  the  sand  and  again 
from  the  water  surface  (E)  are  shown  in  these  figures.  The 
vertical  line  marks  the  beginning  of  the  wave  that  is 
reflected  from  the  top  of  the  plate  and  the  number  gives  the 
arrival  time. 

In  this  case,  the  curve  for  the  arrival  reflected 
from  the  top  of  the  plate  is , 

T  =  0.473  X2  -  32.0  X  +  1020 

and  the  curve  for  the  arrival  reflected  from  the  top  of  the 
sand  is , 

T  «  0.136  X2  -  5.72  X  +  471 

Here,  the  incident  angle  was  28.3  degrees  so  P  was  3.30  *  10  ® 
s/m.  The  ray  parameter  calculations  gave  3.87  cm  for  AX  and 
54  yS  for  dT  so  the  layer  depth  Az  was  3*73  cm.  The  actual 
sand  layer  thickness  was  3.4  cm  so  the  error  was  slightly 
less  than  10%.  The  value  for  sound  speed  in  the  water  by  the 
ray  parameter  method  was  1474  m/s  which  is  about  10%  high. 

3 .  Measurement  3 

In  this  measurement,  an  aluminum  plate  0.1  cm  thick 

was  buried  to  a  depth  of  2.3  ±  0.1  cm.  The  projector  depth  . 

was  8.5  cm  and  the  hydrophone  depth  was  5.1  cm.  The  resuit- 

2  2' 

ing  speed  of  sound  measurements  {determined  by  the  T  -X 
method)  are  summarized  in  fable  III  and  the  T-X  data  and  the 
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2  2 

T *‘~YT  data  are  shown  in  Fig.  4.5.  Figures  4.6  and  4.7 
show  representative  displays  from  the  oscilloscope  at  two 
of  the  ranges  at  which  measurements  were  taken. 


TABLE  III 

Speed  of  Sound  (Measurement  3) 


WAVE 

SPEED  OF  SOUND  (m/s) 

CORRELATION  (%) 

ERROR  (%) 

SAND 

1649 

95.2 

BOTTOM 

1471 

99.9 

1.7 

DIRECT 

1510 

99.8 

4.3 

■ 

The  parabola  for  the 

arrivals  reflecting 

from  the 

top  of  the  plate  was  determined  to  be 

T  =  0.290  X2  -  13.3  X  +  580 

and  the  parabola  for  the  arrivals  directly  reflected  from 
the  top  of  the  sand  is 

T  =  0.0484  X2  +  0.528  X  +  357 

For  the  optimum  geoemtry,  the  incident  angle  was  27  degrees, 
so  that  the  ray  parameter  P  was  3.00  x 10 

The  ray  parameter  analysis  gave  2.53  cm  for  AX  and 
33.38  ms  for  AT;  therefore,  the  sand  layer  depth  was  computed 
to  be  2.42  cm:  less  than  one  percent  from  the  actual  value. 


Figure  4.6.  Measurement  3  at  X  *  21.7  cm. 


Figure  4.7. 


/ 

/ 

The  calculated  value  for  the  sound  speed  in  the  sand  layer 

was  1578  m/s  and  this  differs  by  only  3  percent  from  the 

2  2 

value  obtained  by  the  T  -X  method. 

4 .  Measurement  4 

For  this  measurement,  an  aluminum  plate  0.3  cm  thick 

was  buried  to  a  depth  of  1.8  ±  0.1  cm.  The  projector  depth 

was  8.5  cm  and  the  hydrophone  depth  was  3.5  cm.  The  speed 

of  sound  measurements  is  summarized  in  Table  IV  and  the  T-X 
2  2 

data  and  the  T  -X  data  are  shown  in  Fig.  4.8.1  Representative 
displays  from  the  oscilloscope  at  two  ranges  are  shown  in 
Figures  4.9  and  4.10. 


TABLE  IV 

Speed  of  Sound  (Measurement  4) 


WAVE 

SPEED  OF  SOUND  (m/s) 

CORRELATION  (%) 

ERROR 

SAND 

1582 

99.7 

BOTTOM 

1458 

99.7 

.0.7 

DIRECT 

1461 

99.9 

1.0 

The  parabolic  form  for  the  reflections  from  the  top 
of  the  plate' is 

T  ■  0.0259  X2  +  0.361  X  +  429 

and  the  form  for  the  reflections  from  the  top  of  the  sand  is 
T  *  0.0328  X2  +1.17  X  +380 
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Figure  4.9.  Measurement  4  at  X  *  35.1  cm 


Figure  4.10.  Measurement  4  at  X  *  47.7  cm. 


For  the  optimum  geometry,  the  incident  angle  was  5Q  degrees 

so  that  the  ray  parameter  P  was  5.3  xlG-6. 

The  ray  parameter  analysis  gave  5.80  cm  for  AX  and 

40.4  ys  for,  AT;  therefore,  the  sand  layer  depth  was  computed 

to  be  1.74  cm.  The  calculated  value  for  the  sound  speed  in 

the  sand  layer  was  1647  m/s.  In  both  cases,  the  results 

compare  very  well  with  the  measured  value  for  layer  depth 
2  2 

and  the  T  -X  value  for  sound  speed. 

5.  Measurement  5 

For  this  measurement,  the  projector  depth  was  8.5  cm 

and  the  hydrophone  depth  was  3.5  cm.  The  aluminum  plate 

(0.1  cm  thick)  was  buried  at  a  depth  of  5.0  ±  0.1  cm.  The 

speed  of  souna  measurements  are  summarized  in  Table  V  and 

2  2 

the  T-X  data  and  the  T  -X  data  are  shown  in  Fig.  4.11. 
Samples  of  the  recei  ved  arrivals  at  two  ranges  are  shown  in 
Figures  4.12  and  4.13. 


TABLE  V 

Speed  of  Sound  (Measurement  5) 


WAVE 

SPEED  OF  SOUND  (m/s) 

CORRELATION  (%) 

ERROR 

SAND 

1577 

99.9 

BOTTOM 

1468 

100 

1.4 

DIRECT 

1488 

100 

2.8 

The  parabolic  form  for  the  reflections  from  the  top 
of  the  plate  is 


45, 


and  the  form  for  the  reflections  from  the  top  of  the  sand 
is 

T  =  0,0463  X2  +  0.562  X  +  339 

The  incident  angle  for  determining  the  ray  parameter  was 

33.4  degrees  so  that  the  ray  parameter  P  was  3.70  *10 

The  ray  parameter  analysis  gave  9.59  cm  for  AX  and 

71.5  ys  for  AT;  therefore,  the  sand  layer  depth  was  com¬ 
puted  to  be  4.58  cm — about  six  percent  below  the  actual 
value.  The  calculated  value  for  the  sound  speed  in  the  sand 
layer  was  1904  m/s  which  was  significantly  higher  than  any 
other  of  the  determinations. 

B.  OVERALL  RESULTS 

In  all,  seven  experiments  were  performed  but  the  data 
from  two  of  the  experiments  was  not  good  enough  to  analyze. 
The  other  five  experiments  have  been  reported  above  and, 
except  for  number  5,  the  results  were  quite  Consistent  with 
reasonable  errors. 

Table  VI  gives  a  summary  of  the  results  of  the  five 

experiments  that  were  analyzed  and  compares  the  results  of 
2  2 

the  T  -X  method  and  the  ray  parameter  method.  The  average 
error  in  determining  the  depth  of  the  layer  was  3.5  percent 
so  this  measurement  can  be  done  effectively;  however,  the 
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TABLE  VI 

Summary  of  Brief  Results 
MEASUREMENT 

CONTENT  1  2  3  4  5 


Layer  Depth 
(cm) 

ACTUAL 

EXPERIMENTAL 

ERROR  (%) 

3. 5 (+.05) 

3.56 

0.3 

3.3(+.l) 

3.73 

9.7 

2.3  (+.1) 

2.42 

0.8 

1.8(+.l) 

1.74 

0.0 

5.0  (+ 

4.58 

6.5 

1 

T2-X2  METHOD 

1662 

1647 

1649 

1582 

1577 

4-t  ^ 

RAY  PARAMETER 

1791 

1474 

1578 

1647 

1904 

0  m 

1“ 

DIFFERENCE 

163 

154 

49.6 

18.9  2 

75.9 

8. 

w 

ERROR  (%) 

.  10.0 

9.5 

3.0 

1.2 

16.9 

error  in  the  sound  speed  calculation  (for  the  sand  layer) 

was  fairly  large  (about  8  percent).  Consequently,  the 

sound  speed  determination  is  not  quite  as  reliable. 

In  Table  VI,  the  difference  means  the  values  between  the 

2  2 

average  of  the  five  T  -X  results  and  the  ray  parameter 
method's  value. 
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V.  DISCUSSIONS  AND  CONCLUSIONS 

The  objective  of  this  research  was  to  localize  buried 
objects  (e.g.,  thin  aluminum  plates)  in  water-saturated  sand. 
In  order  to  accomplish  this  objective,  we  needed  to  predict 
the  sound  speed  and  the  depth  of  the  overlying  layer.  The 
results  are  summarized  in  Table  VII. 

2  2 

The  average  speed  of  sound  by  the  T  -X  method  is  1472 
m/sec  for  the  direct  wave,  1446  m/sec  for  the  bottom  reflected 
wave,  1623  m/sec  for  the  wave  through  the  sand. 

Comparirig  the  above  sound  speeds  with  the  value  given  by 
Reference  3  (1481  m/sec  at  20  °C) ,  we  had  0.6%  error  in  the 
direct  wave  speed  and  2.3%  error  in  the  bottom  reflected 
wave  speed. 

Reference  3  gives  the  sound  speeds  1540  m/sec  for  coarse 
silt  and  1730  m/sec  for  quartz  sand.  The  sand  used  in  this 
experiment  is  #30  mesh  sand  (washed,  kilndried  and  graded) 
and  the  measured  sound  speed  is  1623  m/sec. 

We  used  linear  regression  for  finding  the  speed  of  sound 
2  2 

with  the  T  -X  method.  Since  most  correlation  values  of  the 

2  2 

linear  regression  are  over  95%,  the  T  -X  method  is  reasonable 
for  determining  the  speed  of  sound. 

By  using  the  ray  parameter  method,  we  determined  AX,  AT 
and  Az.  The  average  error  between  the  actual  layer's  depth 
and  the  computed  depth  is  3.5%.  The  error  range  is  between 


TABLE  VII 


The  Detailed  Results  of  Total  Experiment 


MEASUREMENT 

1 

2 

J 

4 

5 

DETAHS 

ACTUAL  LAYER 
DEPTH  (cm)  J 

•5(±.05) 

3.3{±.l)  2, 

,3(±.l) 

1.8 (±.l) 

5.0 (±.l) 

■ 

98.3 

94.3 

95.2 

99.7 

99.9 

SPEED  S 

OF 

1661 

1647 

1649 

1582 

1577 

£ 

SOUND 

99.4 

99.6 

99.9 

99.7 

ICO. 

£ 

<N 

(nv/sec)  B 

1416 

1419 

1471 

1458 

1468 

X 

i 

& 

CN 

E-* 

LINEAR 

99.8 

59.9 

99.8 

99.9 

100 

<U 

D 

£ 

COREELA- 

1463 

1437 

1510 

1461 

1488 

TICN  (%) 

1 

i 

P  (10  ^  Sec/m) 

3.32 

3.69 

5.30 

3.70 

AX  (an) 

5.47 

3.87 

2.53 

5.8 

9.59 

£ 

AT  (ysec) 

51.4 

54.0 

33.9 

40.4 

71.5 

<d 

■P 

Az  (an) 

3.56 

3.73 

2.42 

1.74 

4.58 

§ 

ERROR  (%)X 

0.3 

9.7 

0.8 

6.5 

0. 

S 

c.  DIFFERENCE2 

1  (m/sec) 

163 

,  154 

50 

19 

276 

#5 

ERROR  (%)3 

10.0 

9.5 

3.0 

1.2 

'  16.9 

# 

Average  error  of  measuring  Az  =  3.5% 

■  ■ 

* 

Average  error  of 

c^  difference  =  8.1% 

■  • 

^Ihis  error  is  the  value  when  canparing  the  experimental  Az  with  actual 
layer's  depth. 


2  2  2 
c.  difference  is  the  difference  between  c.  of  T  -X  and  the  value  of 

the - J - “ - 1 


ray  parameter  method. 


This  error  is  tr.e  value  whan  canparing  the  c,  di 
average  value  of  T~-X2. 


difference  with  the 


■  ' 
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0%  and  9.7%.  The  ray  parameter  method  is  one  of  the  most 

useful  for  finding  an  unknown  layer's  depth.  But,  the  average 

2  2 

error  of  c1  between  the  T  -X  and  the  ray  parameter  methods 
is  relatively  high  (8.1%). 

In  every  measurement  the  direct  wave ' s  speed  of  sound 

is  always  faster  than  that  of  the  bottom  reflected  wave. 

These  values  should  be  the  same  because  both  are  measures  of 

the  sound  speed  in  water.  This  discrepancy  is  probably  a 

result  of  the  difficulty  in  determining  the  true  horizontal 

distance  between  the  acoustic  center  of  the  projector  and 

the  hydrophone  for  different  angles  of  the  projector.  In 

2  2 

some  cases,  the  difference  in  c^  between  the  T  -X  and  the 
ray  parameter  methods  is  over  10%.  This  may  result  from  the 
choice  of  curve-fit  for  the  T-X  data.  We  tried  several  ways 
to  find  the  best  curve  fit,  but  the  computed  results  were 
quite  different  from  the  expected  ones  except  for  those  from 
parabolic  fit. 

We  had  some  limitations  on  the  measurements  because  of 
the  use  of  a  small  tank.  In  particular,  there  were  inter¬ 
fering  reflections.  So,  some  of  the  signals  arrive  at, the 
same  time,  'that  made  the  analysis  of  the  layer  speed  of 
sound  and  depth  very  difficult.  Also,  the  pulse  length  was, 
too  long  to  be  able  to  resolve  the  thin  layer.  Projector 
ringing  limits  the  pulse  length  so  that  it  cannot  be  made 
shorter  than  a  certain  duration,  because  of  these  problems, 
this  experiment  would  be  easier  to  do  in  a  larger  body  of 


water. 


If  we  do  this  experiment  in  the  ocean  instead  of  small 
water  t=>nk,  then  identifying  the  wave  signals  is  much 
easier  because  of  less  multipath  interference.  In  the  ocean 
we  can  try  A3PER  (airgun  sonobuoy  precision  echo  recorder) 
developed  by  the  Hawaii  Institute  of  Geophysids  as  one  of 
the  recommended  methods  [Ref.  1].. 
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